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Abstract 

We calculate the second-order corrections to the atomic energy level shifts in ordinary 
and muonic deuterium due to virtual excitations of the deuteron which are important for 
ongoing and planned precise experiments in these systems. For light atoms a method can 
be used in which the shift is expressed as integrals over the longitudinal and transverse 
inelastic structure functions of the nucleus. We employ the structure functions arising 
from separable NN potentials of the Yamaguchi and Tabakin form which can be evaluated 
analytically. Special emphasis is put on gauge invariance which requires a consistent 
inclusion of interaction currents and seagull terms. The effect of the D-wave component 
of the deuteron is investigated for the leading longitudinal contribution. We also estimate 
the shift for pionic deuterium. 
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1. Recently the isotope shift of the 2S — IS transition in electronic hydrogen and deuterium 
has been measured with a thirty-fold increase in accuracy compared to previous experiments 
fl]] and prospects are good to increase the present accuracy of 37 ppb (or 22 kHz) by another 
order of magnitude. At this level not only various QED corrections and the finite size of 
the deuteron are important but also virtual excitations of the deuteron (the so-called nuclear 
polarization) cannot be neglected anymore. Using a rough square-well model of the deuteron 
and the closure approximation the IS shift due to virtual Coulomb excitation was estimated 
to be about —19 kHz || which is just the present experimental uncertainty. In light and heavy 
muonic atoms the nuclear polarization shift generally limits the accuracy with which nuclear 
sizes can be extracted since a reliable calculation of these corrections requires knowledge of the 
whole nuclear spectrum ||. The deuteron is unique among all nuclei in that this information 
is available: the quantum-mechanical 2-body problem is solvable and realistic NN potentials 
describe bound and scattering states rather well. Therefore, unlike the usual case, nuclear 
polarization corrections are calculable for the deuteron - presumably with an accuracy at the 
percent level. 

It is the purpose of the present note to evaluate this shift within a more realistic model 
for the deuteron, avoiding the use of the closure approximation and including also transverse 
excitations which have been neglected up to now. For the purpose of a planned experiment 
at PSI with muonic hydrogen M which may be extended to deuterium H we also evaluate 
the corresponding shifts in the muonic case. It should be noted that this does not amount 
to a simple rescaling of the electronic results: since nuclear excitations are of the order of a 
few MeV the electron is higly relativistic in the virtual state whereas the muon can be treated 
nonrelativistically to a good approximation. This leads to different weighting of the individual 
excitations and therefore to quite different results. 

2. We will calculate the nuclear polarization shifts in the atomic S-states of the lepton. These 
are more difficult to evaluate than the corresponding shifts in higher orbits where only the 
longest range multipole is of importance and therefore the only nuclear structure information 
needed is the electric dipole polarizability of the nucleus. In contrast, many multipoles con- 
tribute to virtual excitations from atomic S-states since there is an overlap between leptonic 
and nuclear wavefunctions. In light nuclei, however, the relevant scales (Bohr radius vs. nu- 
clear radius) are vastly different so that the S-wave lepton wavefunction can be considered as 
approximately constant over the nuclear volume 0„o(x) ~ av : the lepton then just acts as a 
static source with four-momentum k = (m, 0) . 

We will evaluate the S-wave nuclear polarization shifts along the lines of ref. || , i.e. not 
by a multipole decomposition but by integrating over the inelastic structure functions of the 
nucleus. The diagrams which contribute to nuclear polarization are shown in Fig. 1. Note that 
we need the "seagull" contribution of Fig. 1 (c) for gauge invariance in a nonrelativistic sytem 
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like the nucleus. Under the mentioned simplifications one obtains 



Ae nO = -^^|0 av | 2 Im / ±jL t flu (q,k)D^(q)D^(-q)f pr (q,-q) (1) 
where 

. / , v _ k-qg^ + jk- q)^ + k p (k - q) v 
(k — q) 1 — m z + ie 

is the leptonic tensor, a = e 2 = 1/137.036 the fine-structure constant and m is the lepton mass. 
Furthermore, D^ p (q) denotes the photon propagator and T pT (q, —q) the forward virtual nuclear 
Compton amplitude. To be more precise, the latter is that part of the full Compton amplitude 
in which the nucleus is in an excited intermediate state. It can be expressed in terms of its 
imaginary part, i.e. by the inelastic longitudinal and transverse structure functions 

Sl/t = E + E Q -E N )\< ^ N \O L/T \^ > | 2 . (3) 

Here uo' = u — q 2 /AM is the internal excitation energy, Eq < the ground state energy and 
Ol/t are the operators for longitudinal and transverse excitations repectively. In Coulomb 
gauge one obtains M 



Ae n0 = -8a 2 R n0 \<j)no(0)\ 2 J Q dq | J du 



K L (q,u) S L (q,u) + K T (q,u) S T (q,u) 



+ R S (q))f(q). (4) 



Here R n o is a correction factor for the variation of the leptonic wave function over the nucleus, 
q is the magnitude of the three-momentum transfer and u the energy transfer to the nucleus. 
The kernels Ki^q^uS) are given in the appendix of ref. for fully relativistic kinematics 
of the lepton. Rs(q) is the contribution from the internal seagull Q. Finally, /(g) = (1 + 
g 2 /(0.71GeV 2 )) -2 describes the electromagnetic formfactor of the nucleon. Actually eq. (4) 
not only holds in the Coulomb gauge but is gauge invariant provided q^T pu (q, —q) = 0. This 
in turn requires current conservation and special conditions for the seagull term || ||| which 
will be discussed below. 

3. Unlike ref. where a phenomenological model for the structure functions of 12 C had to be 
used, the deuteron allows for a consistent calculation of these quantities after a nucleon-nucleon 
interaction has been chosen. For simplicity we take a separable potential of the form (M is the 
nucleon mass) 

V(p,j>') = ~g(p)g(jf). (5) 



1 The center-of-mass seagull term is part of the two-photon recoil correction [jj. Note also that the part of 
the Compton amplitude in which the nucleus remains in its ground state has already been included by solving 
the bound state problem for the lepton with a static potential. 



3 



This is not realistic in the modern sense (it lacks the one-pion exchange tail and all other 
complications of the NN force) but it gives a fairly good description of the low-energy NN 
interaction which should be sufficient for the present accuracy of isotope shift experiments. 
Most important for the present purposes it allows for an analytic evaluation of the structure 
functions. For example, the longitudinal structure function is obtained as [pj] 



S L (q,u>)= /tfWo(p-q/2)| 2 5(u' + E - + Im 

J M 71 

Here the functions C(E) and I(E, q) are given by 

9 2 (p) 



/y + E ,g) ' 
1 + A CV + E ) 



(6) 



drp 



C(E) = 
I(E,q) = 
Note that the internal charge operator is 



ME -p 2 -Me 



d 3 pMp-i/2) 



g{p) 



ME -p 2 + ie 



exp(-q- r) 



(7) 
(8) 

(9) 



where the factor 1/2 arises from the transformation to internal coordinates (a "hat" indicates 
an operator in the nuclear Hilbert space). The ground state wave function is simply given by 



MP) = N 



g(p) 



p 2 - ME 



(10) 



Also note that the first term in eq. (6) is the impulse approximation to the structure function 
whereas the last one describes the final-state interaction. Due to the simple form of eq. (5) 
it only acts in states with angular momentum zero which should be a good approximation for 



low-energy processes. For the Yamaguchi choice |TT 

9y(p) 



1 



p 2 + f3 2 



11) 



all integrals can be performed analytically. Details will be given elsewhere. 

The transverse nuclear polarization shift receives different contributions : first, we have the 
contribution from the transverse (with respect to q ) part of the convection current 



J (conv) (q) 



2M 



P , exp(-q- rj 



(12) 



which does not have a final-state interaction term since all excited states necessarily have at 
least angular momentum one. Second, the spin current 



j (spin) ( 



q) 



2M 



% 1 

H P (<t p x q) exp(-q • r) + //„(<x n x q) exp(--q • r) 



(13) 
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involves the magnetic moments and spins of neutron and proton and gives rise to spin-flip 
excitations. Since the procedure of calculating the transverse structure with these currents is 
similar to the longitudinal case the explicit expressions will not be given here. Note that the 
operators (9), (12) and (13) obey current conservation in the form 



<P 



2\p(q)] - q • (j (conv) (q) + J (spin) (q) ) |p >= o 



(14) 



where T is the kinetic energy operator. 



4. A separable potential of the type (5) is nonlocal and equivalent to a momentum-dependent 
potential 

V = — J d 3 r d 3 x g + ^x^j g — -^x^ exp(— ip ■ x/2) |r >< r| exp(— ip ■ x/2) . (15) 

Minimal coupling p p — ► p p — eA(r p ) in the two-body Hamiltonian H = T + V then produces 
a power series in the electromagnetic coupling constant e 



H = H\ A=0 - e j d 3 yj,(y) A"(y) + - J d 3 y d 3 z My) Aj(z) B^y, z) 

The linear terms due to the potential generate an interaction current 

5V 



(16) 



eAJ fe (y) 



A=0 



SA k (y) 

whereas the second-order terms give rise to the interaction seagull 

5 2 V 



(17) 



e 2 AS -(y,z) = 



5 My) 5Aj(z) 



A=0 



(18) 



With the separable potential (15) one obtains the following expression for the matrix element 
of the interaction current 



<p / |AJ(q)|p>=-^ 7 (Vp + VpO £dug(p-u^j g(p' + (l-u)^ 



2M 



which satisfies 



<P' 



F,p(q) -q-AJ(q)|p>=0 



(19) 



(20) 



so that all together our currents are conserved under the time-evolution of the full Hamiltonian. 
Similarly one obtains the following matrix element of the interaction seagull 

< p'| ASy(q) |p >=< P f = 0,p'| J d 3 xexp(iq-x) A%(x,0) |p,P, = > 
A /H-i 



4M 



(V p + V p/ ) t (V p + V p/ ) j du (1 -\ u \)g(p- u^J g (p' - u^j (21) 
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which comes in addition to the kinematical (internal) seagull 5y/2M. 

It is well known || |9[ that for full gauge invariance of the Compton amplitude a relation 
between the interaction current and the interaction seagull is needed. In the present case where 
no overall momentum is transferred to the target it reads 

<p'| [p t (q) ) AJ jfc (q)] |p>=g < p'| AB 3k {q) |p > . (22) 

With the interaction current and interaction seagull as given in eqs. (18) and (21) this re- 
lation can be shown to be fulfilled. Thus the nuclear polarization shift is gauge invariant if 
the interaction terms in current and seagull are taken into account consistently. The seagull 
contribution in eq. (4) is now explicitly given by 

Rsiq) = 1 8^T (I ~ JrJ+q* ) (23) 



K(q) = ^ fdu^ \h 2 (0) - h 2 {uq)] (24) 



where 



q" Jo m z 

h(x) = J d 3 p MP) 9 (p - f ) ■ (25) 

5. We now turn to the numerical results of our calculation. We have evaluated the nuclear 
polarization shift consistently with the Yamaguchi separable form (5) using the value (3 = 286 
MeV. We have checked numerically that our structure function fulfills the non-energy-weighted 
sum rule 

poo 

/ dojS L (q,u) = 1 - F 2 (g) (26) 
Jo 

to better than 1 part in 10 5 . Here Fo(q) is the elastic formfactor calculated directly from 
the ground state wave function (10). The electric dipole polarizability (which is not the only 
relevant quantity for S-wave shifts) was found to be 0.613 fm 3 compared to experimental 
values of 0.61 — 0.70 fm 3 The point rms-radius in this model is 1.92 fm compared to the 



experimental value 1.96 fm. This shows that the simple Yamaguchi parametrization describes 
the deuteron properties and therefore the low-energy triplet NN-interaction reasonably well. 
For the spin-flip excitations one also needs the parameters in the singlet channel. Again 
following Yamaguchi we assume (3 S = fit and determine the corresponding strength parameter 
from the singlet scattering length a s = —23.69 fm. 

For non-relativistic point hydrogen wave functions one has 

l«0)! 2 = 4^ (27) 
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where as = l/(am re d) is the Bohr radius and m re d the reduced mass of the lepton. Writing 



the shift Ae is then independent of the atomic state. We have evaluated the double integral in 
eq. (4) by Gauss-Legendre numerical integration with up to 3 x 72 points. Our results for the 
different contributions and for the total Ae are listed in Table 1. It should be emphasized that 
the integrand from the transverse convection current has a 1/q - divergence for small q which 
, however, is exactly cancelled by the seagull due to the gauge condition for the two-photon 
operator. This can also be seen from eqs. (4) and (23) : at low q the kernel KT{q,uj) behaves 
like to — l/Amqu |J and Siegert's theorem tells us that 



where k = k(0) is the dipole enhancement factor (k = 0.176 for the Yamaguchi potential). 
However, the resulting contribution to the energy shift is exactly opposite in sign to the q — > 
0-limit of the seagull contribution (23). We have checked numerically that our transverse 
structure function fulfills the sum rule (29) to sufficient accuracy. Note that in ref. || the 
seagull was represented by eq. (29) for all q; we have found numerically that the extra in- 
dependence only gives a tiny contribution to the full shift. 

In Table 1 we therefore give only the combined contribution of transverse convection and 
seagull excitations. It is seen that it is bigger in electronic deuterium than in muonic deuterium 
because the electron velocity is higher in the first case. As the spin current contribution vanishes 
for q = (see eq. (13)) it can be given separately. However, numerically it turns out to be of no 
great importance. The same can be said of the interaction terms which nearly cancel the spin 
contribution. The smallness of the interaction terms is welcome since the nonlocality of the 
Yamaguchi separable potential is somehow artificial and only partly simulates exchange current 
effects. It should be kept in mind that the individual contributions are gauge-dependent and 
that only the total Ae is a meaningful physical quantity. The size of the transverse and seagull 
terms , however, indicates the errors one usually makes if only the longitudinal excitations are 
taken into account. As to the numerical accuracy, we have checked that the results in Table 1 
are accurate to one part in the last digit. 

In order to estimate the model dependence of these results we also have calculated the 
dominating longitudinal contribution for the Tabakin separable potential JOJ which describes 
both attraction at low energies and repulsion at higher energies in the S-wave phase shift. 
The principal value integrals in eqs. (7) and (8) have now been performed numerically which 
considerably increased the computing time. Again the sum rule was checked and an electric 
dipole polarizability of 0.623 fm 3 was obtained. The values for the shift given in Table 2 are 
estimated to have an accuracy of about three parts in the last digit. Despite the different 



Ae 



nO — 




n 



(28) 




(29) 
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functional parametrization of g(p) the result for Ae is very close to the Yamaguchi one which 
shows that only low-energy properties of the NN-interaction are important for the nuclear 
polarization shift. 

Finally we also have investigated the influence of the D-state admixture in the deuteron by 



using the Yamaguchi separable potential with tensor force [14 



g(p) = 9y(p) + -4 t y(p) S Pn(p) (30) 



with 



Trip) = -tJ^. (31) 



(p 2 + 7 

S pn (p) — 3— — ^ 2 n - — er p -cr n . (32) 



The original parameter values of ref. [14j lead to an asymptotic D/S ratio of 0.0285 which is 



quite reasonable when compared with modern values [13] . The dipole polarizability is calculated 
to be 0.625 fm 3 . Since the algebra including the tensor force is more involved the sum rule 
check (to one part in 10 3 ) is nontrivial. 

To convert these numbers to the actual nuclear polarization shifts for the S-levels we need 



the finite size correction factors R n o . Using the approximate atomic wavefunctions of ref. |16 
one obtains in first order in the ratio nuclear radius/Bohr radius 

< r 2 > l/2 

R n0 ~ 1 - 3.06 . (33) 

a B 

The numerical factor in this equation was determined by evaluating the ratio of various mo- 
ments of the charge distribution with the Yamaguchi ground state wavefunction. Eq. (33) 
gives n-independent correction factors = 0.9793 and = 0.99989. Of course, on the 
present level of accuracy one can practically neglect these correction factors. 

We estimate the accuracy of our theoretical predictions in the following (rather conservative) 
way: the accuracy of the calculated longitudinal shift is taken as three times the model- 
dependence shown in Table 2 and we assign a 20 % error to the transverse current contribution 
and 50 % one to the interaction pieces. Adding these errors linearly we therefore arrive at the 
final result for the nuclear polarization shifts in electronic and muonic deuterium 

AeS = (-20.5 ± 1.3) — kHz (34) 

n 6 

Ae^ = (-11.6 ±0.5) — meV. (35) 

n 6 

We also have estimated the nuclear polarization shift in pionic deuterium by replacing the muon 
mass by the pion mass. The longitudinal and the transverse convection current contribution of 
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the present formalism should give a reasonable value also for a heavy spin zero particle because 
to a good approximation it can be treated nonrelativistically with no difference between a Dirac 
and a Klein - Gordon description. In this way we obtain 



AeS ~ -28 A, meV . (36) 



1 



n 



For n = 1 this is a factor of two smaller than the precision aimed at in an ongoing experiment 
at PSI to measure the strong interaction shifts in pionic hydrogen [HJ . 

If the future isotope shift experiments in electronic deuterium actually reduce the experi- 
mental accuracy to about 1 kHz it would be worthwhile to repeat the present calculation 
with a modern NN-potential like the Paris potential ||18|| . Also it may be useful to employ 
a general gauge like the "a-Lorentz gauge" [[0| which interpolates between Coulomb and 
Lorentz gauge thereby demonstrating explicitly the independence of the numerical results from 
the gauge parameter a. Finally, before discrepancies between theory and experiment in the 
isotope shifts are taken serious one should include second-order effects also in the analysis 
of electron-deuteron scattering experiments which extract the root-mean square radius of the 
deuteron. In the case of 12 C similar discrepancies between electron scattering data and muonic 
energy shifts seem to disappear [OT when second-order effects are taken into account in the 
analysis of both experiments. 
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Table 1 : 

Contributions to the nuclear polarization shift Ae ( see eq. (28) ) for electronic (e) and 
muonic (/x) deuterium in Coulomb gauge. The Yamaguchi S-wave separable potential has been 
used throughout. The different contributions are labeled by L : longitudinal; T ( - conv ' ) + S : 
tranverse convection current + seag ull; T( spin ) : transverse spin current; A(T + S) : interaction 
transverse current + interaction seagull. 



Contribution 


e [kHz] 




meV] 


L 


-18.31 




-11.77 


rp(conv) _j_ g 


- 2.25 




- 0.06 


rp(spin) 


+ 0.33 


+ 0.03 


A(T + S) 


- 0.31 




-0.02 


total 


-20.54 




-11.82 



Table 2 : 

Longitudinal nuclear polarization shift Ae for different separable NN-potentials. 



NN-Potential 


e [kHz] 


H [meV] 


Yamaguchi S-wave 
Tabakin S-wave 13 
Yamaguchi S+D-wa 


11 1 

ve 


14| 


-18.31 
-18.54 
-18.45 


-11.77 
-11.92 
-11.86 



Figure caption 
Fig. 1 : 

Second-order contributions to the nuclear polarization shift: (a) box graph, (b) crossed graph, 
(c) seagull graph. 
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